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^ \ Abstract: We investigate in detail the structure of mesonic vacua of N=l U{N C ) super- 

symmetric gauge theory with Nf flavors from the matrix model. We show that the Witten 
index from the matrix model calculation agrees with a result from field theoretical analy- 
sis. We also discuss the relationship between a diagrammatic summation and direct matrix 
integration with insertion of a variable changing delta function. Using this formalism, we 
obtain the quantum moduli space and evidence of the Seiberg duality from the matrix 
^ , models. 
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1. Introduction 

Dijkgraaf and Vafa recently showed that an effective superpotential of N=l supersymmetric 
gauge theory is obtained from a hermitian matrix integral over a tree level superpotential 
EL @» Ml ■ This observation is originally formulated by using a string theoretical realization 
of gauge theory, but even apart from string theory the Dijkgraaf- Vafa (DV) procedure 
itself produces field theoretically a powerful principle to calculate the effective glueball 
superpotential. After their works, a lot of papers, which discuss from various point of 
view, have rapidly appeared in P|- p^j . 

The model discussed in ||, ^, pj is mainly the N=l gauge theory coupled to adjoint 
chiral superfields with a general superpotential of a polynomial. And also N=l* and Leigh- 
Strassler deformed theory from N=4 theory are discussed in detail by [j| EM. On the 
other hand, theory with chiral fields in a fundamental representation, which is the so-called 
supersymmetric QCD (SQCD) theory, has much interesting dynamics like a generation of 
the Affleck-Dine-Seiberg (ADS) superpotential Seiberg duality (3(], [32|, etc. So we 
expect that the matrix model technique sheds new light on the non-perturbative dynamics 
of N=l SQCD and phenomenologically realistic models. 

A breakthrough in the calculation of matrix models with flavors has been achieved by 



12]. They find that for U(N C ) gauge theory with Nf flavors in the fundamental repre- 
sentation, the matrix integral is completely performed and infinite series in a perturbative 
expansion could be represented by a single analytic function. Using their results, a genera- 
tion of ADS superpotential is discussed in [14, 15|. However, these matrix model calculation 
has a discrepancy in a sense of the number of supersymmetric vacua (Witten index) com- 



pared with results of [^] except for the case of Nf = 1 [14|. This implies that we overlook 
some information in the matrix model calculation. 



Moreover, the generation of the ADS superpotential is discussed in a different way [16] 
by insertion of a matrix valued delta function, which is called as the Wishart integral in 
the mathematical literature |34], |3{|. However, the relationship between these two different 
methods is still unclear. If we naively use the result from the Wishart integral for the 
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N c < Nf case, the matrix integral indicates an existence of the Seiberg duality 21]. 
But we need more careful analysis for this attractive problem. 

The purpose of this paper is to resolve the discrepancy of the Witten index between 
the effective mesonic field theory and matrix model calculation, and to find an exact cor- 
respondence between the perturbative series expansion and the direct matrix integration 
with insertion of the delta function. The former is just done by perceiving ambiguity of 
the functional representation in a perturbative infinite series. The latter can be found by 



refining the arguments in [16] from the first principle of the matrix integral. In our argu- 
ment, we obtain the exact ADS superpotential without the shift of the mass term for the 
meson field. 

The paper is organized as follows. In the subsequent section, we extend the analysis 
of the effective mesonic vacua which is discussed in |33| for the SU(N C ) case. This analysis 
tells us a disagreement of the Witten index with the previous matrix model calculation. 
We reproduce the matrix integration for a general mass matrix in section 3 and we show 
one to one correspondence between the vacua of the meson and glueball effective theory. 
In section 4, we discuss in detail the relationship between the results of the perturbative 
sum of the matrix integration and one from the direct matrix integration. We will give a 
meaning of the mesonic effective theory in terms of the matrix models. 

2. Vacuum structure of effective meson theory 

We consider U (N c ) gauge theory with Nf flavors in the fundamental representation in the 
following. The tree level superpotential of this theory is obtained by a flow from N=2 
SQCD theory with fundamental flavors Qf and Qf by adding a mass \x for the adjoint 
scalar in N=2 vector multiplet 

W tiee ($, Q, Q) = »Tr$ 2 + V2Qf$ b a Qi + rmQ^Ql (2.1) 

where i, j = 1, . . . , Nf stand for the flavor indices and a,b = 1, . . . ,N C stand for the color 
indices. 

If we integrate out $ first, we expect the effective superpotential for mesons from the 
principle of holomorphy and dimensional counting as follows 

Weff(M) = "^TrM 2 + Tr(mM) + Wads> (2.2) 
where we define the gauge singlet meson variable by M\ = QIQa and IVads is an instanton 



induced effective ADS superpotential proposed by [29| 



/ .Wc-Nf \ Nc-N f 

w ^=^- N '\-kM-) • (23) 

with the scale of the N=l SQCD theory Ai. The vacuum structure of this effective theory 
is carefully analyzed in f3lf| for SU(N C ) theory. We follow their arguments below for the 
U(N C ) gauge group. 
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Using the U(Nf) global symmetry, we can diagonalize M as 

M = diag(Ai,A 2 ,...,A JV/ ). (2.4) 



Substituting the diagonalized M into the effective superpotential (2J2), we have the vacuum 
equations for each Aj 

x A {3N c ~N f )/(N c -N f ) 

+ mi ~ (n^w-*/) A - rl = °' (2 - 5) 

and equivalently this reduces to a quadratic equations 

Xf - nmiXi + fiX = 0, (2.6) 

where we define 

X = A f^/)/(^-^/) ( JJ A .)l/(iV / -Ar c )_ (2J) 

5 

These quadratic equations have two different solutions for each \. We now introduce 
two indices i and i' which run for solutions with negative and positive signs, respectively, 
within Nf indices (this choice is for later convenience), that is, we choose 



\ { = - - Aa { Xj , i runs for r within N f (2.8) 

Aj/ = (l + y/1 - AavX} , i' runs for N f -r within Nf, (2.9) 

where r is the number of the solutions with negative signature and we define a,- = — 

We now recall that X is defined in in terms of Aj's. So putting ( |2.8j ) and ( fO| ) 
into the definition of X, we obtain the following equation for each fixed r 



3N C -N f 



X = A*-"' (det (^)) [J (l - yi^x) ^ [J (l + xA - 4a,X 

i i' 

(2.10) 

and reduce to 

x»- = ^jjA 3Nc n (i + v/r^x) n (i - yr^ia^) , (2.11) 

i i' 

where A 3Nc = (det m)A? c Nf is a dynamical scale of the pure U(N C ) gauge theory. 

We consider in the following the case that the mass matrix is proportional to an 
identity in order to see the explicit structure of solutions in eq. ( |2.11 ). In this simple case, 
the general solution of the eigenvalues preserves a global U(r) x U(Nf — r) symmetry and 
eq. fl2.11p becomes 

X N ° = ^jA 3N < (l + Vl - 4aA:) r (l - Vl - 4«x) 7V/ ~ r , (2.12) 

where all ctj's are identical to a now. If we assume r > Nf—r without loss of generality, 
we can rewrite the vacuum equation into 

Y N °~ r (l - Y) N °~ N f+ r = A SN °a N «, (2.13) 
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where we define Y = (1 + \/l — AaX)/2. 

We can see that eq. (gig ) is a polynomial of degree 2N C —Nf in y. So we expect that 



eq. ( p.12 ) has 2N C —Nj solutions for each fixed r as well as eq. ( 2. 13| ) . The number of ways 



for choosing r negative signature is given by N f C r , but if we note that the polynomial 
( 2 . 1 3[ ) is invariant under exchanging the signs of r positive roots with Nf—r negative ones 



in eq. ( I2.12I ) 1 , the essential number of choice is a half of these. So we can count, in terms 
of X, the total number of vacua (Witten index) as 



N f 

y /' > V.. - Yrl^ 

(2N C - N f )2"f~ x . (2.14) 



iV f 

Tr(-1) F = h 2 N c -N f )J2 Nf C r 



r=0 



In comparison with the results in |15| , there exists a discrepancy of the Witten 
index with the factor of 2 N f~ 1 except for the Nf=l case [14]. We will refine the matrix 



model calculation in the next section to resolve this discrepancy. 
3. Matrix integral over flavors 

We will discuss the derivation of the glueball effective theory from the matrix model pro- 



cedure given by 0, H]. We here extend the diagrammatic calculation with flavors in [12] 
to a general mass matrix. 



As conjectured in [12], only the diagram with one flavor index loop contributes to 
the effective superpotential. A topology of this diagram is a disk with one boundary and 
we have to sum up all other color index holes. We can insert S = TrA 2 into each color 
index holes except for one index loop in order to cancel fermionic zero modes. In the 
diagram which contains only the color index loop, we can choose one hole in which no S 
operator is inserted among total h holes. Therefore, the contribution to the superpotential 
is proportional to N c hS h ~ l , see Fig. 1, where h in the prefactor is the number of choices 
of no operator inserted loop and the factor N c comes from the trace of Chan-Paton color 
index. This means that the contribution to the effective superpotential from the free energy 
of the matrix model is N c ^^^-, where J- X= 2{S) is a sum of the planar diagram with the 
Euler number x = 2, or the genus zero. 

On the other hand, if we consider the diagrams with one flavor index loop (boundary), 
we can not insert the operator S into the flavor loop. So we must insert into h — 1 color 
indexed holes and there is no other choice of the insertion. This gives just the contribution 
of S h ' to the superpotential, where h! = h — 1. So the contribution to the superpotential 
directly comes from a sum of the diagram with one boundary ■7 r x= i(S') instead of the 
derivative of the free energy. We finally take a trace over the flavor index. This leads 
to the factor Nf if the system possesses a U(Nf) global symmetry, but it becomes a sum 
over the different function with each flavor indices in general, namely Yli=i -F^iiS), where 
each function depends on some parameters, which break the global symmetry, like different 
masses. 



1 There is no symmetry like this for the general mass matrix. 
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N c hS h - 1 N f S h ' 



Figure 1: Open string world sheet diagrams of the flavor theory that contribute to the superpo- 
tential. We denote color index loops as solid lines and a flavor index loop as a dotted line. There is 
no operator S insertion in the outer index loop. The contribution from the diagram with one flavor 
index loop (one boundary) is proportional to NfS h . 



As mentioned above, if we use a mass matrix which is proportional to an identity 
matrix and preserves an U(Nf) global symmetry, we have only the contribution from the 
diagrams with one boundary is always proportional to Nf. This hides the detailed vacuum 



structure from the matrix model calculation. So we extend the calculation of [12] to the 
general mass matrix, which is diagonalized by the global symmetry at least, to compare 
with the field theoretical results in the previous section. 

Now let us consider the matrix integral with the tree level superpotential (2.1) to 
evaluate the contribution from the diagrams with one boundary, 

1 fr^r^r.A, f 



V6l(U(N c )) 



J [d$}[dQ][dQ]exp^-^Wtrec($,Q,Q)} , (3.1) 



where Yol(U(N c )) is a volume of the gauge group U(N C ). We regard the quadratic terms 
of the matrices in Wt Tee (&, Q, Q) as a kinetic term and the Yukawa interaction term as 
an interaction of the perturbative expansion. We find two point correlators (propagators) 
from the tree level superpotential as 

The contribution to J- X= 2{S) comes only from the volume factor and the measure of 
a gaussian integral of $ since W tree does not include any self- interacting part of <I>. So the 
matrix integral fl3.1|) has the following series expansion with ignoring the volume factor 



z = P o Wy. (^) ((0*0)1(0*0)3 • • • (0*0)»>- ( 3 - 3 ) 

Using the Feynman rules obtained from the superpotential, we can find weights of the 
connected diagrams with one boundary and contract each two point correlators, thus we 
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get 



^ (2fc-l)! (V2N C \ 




Exchanging an order of the summation, we finally obtain 



(3.5) 



where on is the same one defined in the previous section. 

We now consider the series in (3.5) which is written by a generalized hypergeometric 
function and has a branch cut discontinuity running from | to oo in the complex plane. 
This hypergeometric series could be represented by a standard function as 



y (2fe-l)! fc _l 



i 



2/(±)(x) 



log /(±)(x), 



(3.6) 



where f^'(x) are two solutions of the quadratic equation f 2 —f+x = 0, namely, defined 
by 

/(±)(i) = l±vT^ p ?) 

The sign of f^\x) depends on which branch x exists. 

If we use the above functional representation of the series ([O]), we obtain the free 
energy from the connected diagrams with one boundary 



F x =i{S, ati 



N f , 



log f&XatS) 



(3.8) 



2/(±)( ai 5) 

where the signs are suitably chosen for each i. We notice that the vacuum of [12, 15 
represents only the case that all signs are the same, namely the r=0 case. This structure 
of vacua is preserved even though the mass matrix is proportional to the identity matrix. 

Therefore finally we obtain the effective superpotential from the above matrix integra- 
tion in terms of the glueball field 



W cS (S) = N c d ^ iS) 



dS 



S 



N r . 



i-log^-gQ-^-^-log/WKS)) 



(3.9) 



The F-flatness equation with respect to S gives 

s»* = ^A^n i 1 + vt^s) n (i - v^^s 



(3.10) 
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where i and i' run for index of the positive and negative roots, respectively. 

This equation is exactly equivalent to eq. ( 2.1 1| ) identifying X with S. So these theories 



have obviously the same structure of the vacuum solutions. If we especially consider the 
case that all c^'s are identical, the vacuum equation has 2N C —Nf solutions with each fixed 
choice of the signs and the number of the independent way to choose the signs is 2 N f'~ 1 for 
the same reason mentioned in the previous section. So we conclude that the Witten index 
of this theory is 

Tr(-1) F = (2N C - iV / )2 Ar /~ 1 , (3.11) 



which agrees with ( 2.14 ). 



We note here that the above vacuum structure still holds in the limit of m; — ► and 
fi — > oo with CKj fixed, namely we can take a massless limit of quarks without the adjoint 
scalar fields <3? in N = 1 theory. However if take the limit of \i — > oo first, which corresponds 
to ai — ► 0, then the vacuum equation reduces to (no more choice of the signs) 

S N c = (detm)Al Nc - Nf , (3.12) 

which gives only N c vacua as well as the pure U(N C ) case. We must be careful that these 
limits are not commutable. 



4. Complete derivation of the ADS superpotential from matrix model 

In this section, we will see more explicit correspondence between the effective meson theory 
and the direct matrix model calculation with a complete proof. We again begin with the 
matrix integral 

Z = Yol ^ Q m) J [d*][dQ][dQ]exp|-lw' tr ee(*,Q,Q)| , (4.1) 

where Mq = (det(m/27rg s )) Nc is a normalization constant which is determined by the 
gaussian integral over Q and Q. 

To compare with the effective theory written by the mesonic variables M? = QfQa, we 
consider a change of variables in the superpotential by insertion of a matrix valued delta 
function 

e _^ trcc($ , Q ,Q) = / [dM]6{M _ gg )e -i^(*,M) ] (42) 

where 

W tree ($, M) = ^Tr4> 2 - — TrM 2 + Tr(mM), (4.3) 
with <|k = & b a + -Ji^QiQa- The matrix valued delta function is defined by 

5(M - QQ) = — / [dT] exp \iTrT(M - QQ)\ , (4.4) 

where T is an Nf x Nf matrix. Note that this variable change is valid for N c > Nt since 
each matrices have the same degrees of freedom. If N c < Nf, then we have to introduce 
more constraints of variables like baryons. 
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Substituting (|4.2|) into ( |4. l|) and integrating gaussian integrals of Q, Q and <I>, we find 

-5—3 / [d$][dQ][dQ][dA/][dT] exp (-— WW*, Af) + iTr T(M - QQ) 

(2vr) iV /Vol(C/(iV c )) V I 9s 

^——7 / WW (det(iT)T Nc exp {--Wtree(M) + *Tr Tm) , 



(4.5) 



where Wtree(Af) = Tr ^— j^M 2 + mM^j is the tree superpotential of the mesonic part only. 
Let us consider the following integration over T first. We define an integral 

I Nc , Nf (M) = J [dT] (det(iT))-^ e 4T ™. (4.6) 

To perform this integral, we decompose the matrices T and M into smaller parts 

*-(;?)• -(:£)• 

where r and /0 are scalars, t and t> are iVy — 1 vectors and T' and M' are {Nf— 1) x (Nt— 1) 
matrices. This decomposition leads to the decomposition of the determinant of iT 



det(iT) = Ur + t^iT') -1 *) det(iT'). (4.8) 
Now we can integrate over the real scalar r by a residue integral 

I Nc ,N f {M) = ^- p ^J[dT'] (det(z'T')) ~ Nc exp {iTr T'M'} 

x y dtdt t exp|- / 9(tt + at)(iT / )" 1 (t + a)-i^( iT ') l ,| i (4.9) 

where a = —^(iT')v. Next if we notice that the integration over t and r i s a gaussian 
integral with a constant shift, we obtain the following recursive relation 

lN c ,N f (M) = (2 "^V ^ / [dT'] (det(,T'))^- 1 exp j.TrT' (V - I„„t) } 

= ^r ^ 1 ^- 1 ( m/ ■ ^ vt ) • (4io) 

We can easily solve the above recursive relation from the initial condition Jjv c _jv>+i,i(p) - 
v(Nc % f+1) P Nc ~ Nf if and only if N c > N f . So we find 

^ W = zzWt , : ( det M ) c f 

n&w-j + i) 

= (2vr)^/ + ^ G2(Ar r iV/ | 1) (det M) N °~ Nf , (4.11) 
v ; G 2 (-/V c + 1) v ; v ; 
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where G 2 (z) is the Barnes function defined by 



G 2 (z + l) = T(z)G 2 (z), G 2 (l) = l, (4.12) 

and which has an asymptotic expansion at the large z 

z 2 ( 3 \ z 1 

log G 2 (z + 1) ~ — (log z - -J + - log 2tt - — log z. (4.13) 

We would like to emphasize here that the formula ( 4.11] ) does not contain a term 
e -TrM j n con t ras t w ith the discussion in flq , 34, pfifl . This term has a difficulty in a field 
theoretical sense since it is proportional to a mass term of the meson field and shifts the 
meson mass. However, there is no problem in our arguments. 

Now if we recall the volume of the U(N) gauge group is given by 36, 37] 

(2ir)3 A " + i A ' 

Vol(t,(A '» = -GWTW <" 4) 

plugging back the formula ( |4.11| ) into eq. ( f4.5|) , we obtain 



2 ~ VWM-N,)) j l " M] (d6t M)N " N ' e ~" W '- m - (415) 

The volume factor is changed from U(N C ) to U(N C — Nf). This means that the Veneziano- 
Yankielowicz term is modified to (N c — Nf)S log S due to the matter contribution. 

We first consider to complete the above integration over M. In the large N c , Nf limit 
with Nt/N c fixed, the above integration over M can be evaluated by the saddle point 
approximation. If we again diagonalize M like ( |2.4| ) 2 , the integral of M is dominated at 
the solution of the saddle point equations 

--Xi + mi- SX; 1 = 0. (4.16) 

This quadratic equations has two solutions for each Aj 

(±) f l±Vl-4a l S \ 
A- ' = iimi[ . (4.17) 



Picking terms proportional to l/g s from the approximated M integral part and the volume 
factor in — log Z, we find exactly 

N f / 

^i(S,ai) = -SJ2 [I ~ 2fi ±) {aiS) ~ lQ g / (±) («^)J . (4-18) 

where f^(x) are defined in ( |3.7| ) and signs are properly chosen for each i. This final 
result, of course, is the same as ( |3.8| ) in the previous section since we have just performed 
the matrix integration over all variables in the alternative way. 

2 This does not bring about the Vandermonde determinant since M is always diagonalized by the fixed 
U(Nf) global symmetry. 
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We have here integrated all matrix variables to the end, so we have skipped the in- 
termediate mesonic effective theory with the ADS superpotential. To obtain the ADS 
superpotential from this matrix integral, we have to forget temporally the integral over M 



in (4.15) and use the DV procedure before the integral of M. We can justify this procedure 
since the size of M depends only on Nf and the M integral does not contribute to the 
function of S. We find the expansion of the free energy from — log Z and get 



W cS (S, M) = N c 8:Fx -f S) - 2mr Q S + F x =i(S) 

= S [N c (l-log(S/A 3 )) - N f (1-logS) - log (det(mM))] + W tree (M). 

(4.19) 



Next we consider the minimaizing of S, namely if we solve the F-flatness condition with 
respect to S and substitute into the solution to ( f4.19| ), we obtain the effective meson theory 

W eS (M) = W tiee (M) + W ADS (M), (4.20) 

where Wads is exactly the same as (|2.3| ) with the identification of the scale A± Nc Nf = 



A 37Vc /detm. If we minimize the superpotential fl4.20| ), we apparently see the same vacuum 



structure as the case where all matrix variables are integrated since we just minimize with 
respect to M and S in each case. 

Now we consider the special case of N c =Nf. We are able to use the above argument 
for N c =Nf excluding the baryons, B = detQ an d B = detQ. All logarithmic functions of 
S disappear from eq. ( |4.19| ) in this case. Then we have 

W eS (S, M) = Slog I J + W tree . (4.21) 

From the F-flatness equation of S, we obtain the following equation without any gluino 
condensation 

detM = A 2 1 Nc . (4.22) 

Therefore we obtain a quantum modification of the classical moduli space from the matrix 
model H when B = or B = 0. 

Finally we note that the integration of the delta function exchanges the volume factor 
of U(N C ) for one of U(N C —Nt). If we naively apply the formula ( 4.1 1| ) to the case of 
N c < Nf, we can rewrite eq. ( |4.15| ) up to some constants as 

Z ~ ——I Ar „ / [dM][dq][dq]e-^ iWt ^ {M)+ " Mq) , (4.23) 



Vo\(U(N f -N c )) 

where q and q are (Nf—N c ) x Nf matrices, namely, which represent the fundamental flavors 
of U(Nf—N c ) gauge theory. This fact strongly suggests the appearance of the Seiberg 



duality from the matrix models [20, 21|. Moreover, we comment that the above rough 
derivation of Seiberg duality is valid even without any mass deformation. However, the 
degrees of freedom of the quark matrices Q and Q can not be replaced by only the mesonic 
variables for N c < Nf. So we should need more careful analysis with some constraints or 
baryon matrices |26[|. 
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